\(\int x^4 (a+b \sec ^{-1}(c x)) \, dx\) [3]

Optimal result

Integrand size = 12, antiderivative size = 89 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=-\frac {3 b \sqrt {1-\frac {1}{c^2 x^2}} x^2}{40 c^3}-\frac {b \sqrt {1-\frac {1}{c^2 x^2}} x^4}{20 c}+\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {3 b \text {arctanh}\left (\sqrt {1-\frac {1}{c^2 x^2}}\right )}{40 c^5} \]

[Out]

1/5*x^5*(a+b*arcsec(c*x))-3/40*b*arctanh((1-1/c^2/x^2)^(1/2))/c^5-3/40*b*x^2*(1-1/c^2/x^2)^(1/2)/c^3-1/20*b*x^
4*(1-1/c^2/x^2)^(1/2)/c

Rubi [A] (verified)

Time = 0.03 (sec) , antiderivative size = 89, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.417, Rules used = {5328, 272, 44, 65, 214} \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {3 b \text {arctanh}\left (\sqrt {1-\frac {1}{c^2 x^2}}\right )}{40 c^5}-\frac {b x^4 \sqrt {1-\frac {1}{c^2 x^2}}}{20 c}-\frac {3 b x^2 \sqrt {1-\frac {1}{c^2 x^2}}}{40 c^3} \]

[In]

Int[x^4*(a + b*ArcSec[c*x]),x]

[Out]

(-3*b*Sqrt[1 - 1/(c^2*x^2)]*x^2)/(40*c^3) - (b*Sqrt[1 - 1/(c^2*x^2)]*x^4)/(20*c) + (x^5*(a + b*ArcSec[c*x]))/5
- (3*b*ArcTanh[Sqrt[1 - 1/(c^2*x^2)]])/(40*c^5)

Rule 44

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[(a + b*x)^(m + 1)*((c + d*x)^(n + 1
)/((b*c - a*d)*(m + 1))), x] - Dist[d*((m + n + 2)/((b*c - a*d)*(m + 1))), Int[(a + b*x)^(m + 1)*(c + d*x)^n,
x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && ILtQ[m, -1] &&  !IntegerQ[n] && LtQ[n, 0]

Rule 65

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 272

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a
+ b*x)^p, x], x, x^n], x] /; FreeQ[{a, b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]

Rule 5328

Int[((a_.) + ArcSec[(c_.)*(x_)]*(b_.))*((d_.)*(x_))^(m_.), x_Symbol] :> Simp[(d*x)^(m + 1)*((a + b*ArcSec[c*x]
)/(d*(m + 1))), x] - Dist[b*(d/(c*(m + 1))), Int[(d*x)^(m - 1)/Sqrt[1 - 1/(c^2*x^2)], x], x] /; FreeQ[{a, b, c
, d, m}, x] && NeQ[m, -1]

Rubi steps \begin{align*} \text {integral}& = \frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {b \int \frac {x^3}{\sqrt {1-\frac {1}{c^2 x^2}}} \, dx}{5 c} \\ & = \frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )+\frac {b \text {Subst}\left (\int \frac {1}{x^3 \sqrt {1-\frac {x}{c^2}}} \, dx,x,\frac {1}{x^2}\right )}{10 c} \\ & = -\frac {b \sqrt {1-\frac {1}{c^2 x^2}} x^4}{20 c}+\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )+\frac {(3 b) \text {Subst}\left (\int \frac {1}{x^2 \sqrt {1-\frac {x}{c^2}}} \, dx,x,\frac {1}{x^2}\right )}{40 c^3} \\ & = -\frac {3 b \sqrt {1-\frac {1}{c^2 x^2}} x^2}{40 c^3}-\frac {b \sqrt {1-\frac {1}{c^2 x^2}} x^4}{20 c}+\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )+\frac {(3 b) \text {Subst}\left (\int \frac {1}{x \sqrt {1-\frac {x}{c^2}}} \, dx,x,\frac {1}{x^2}\right )}{80 c^5} \\ & = -\frac {3 b \sqrt {1-\frac {1}{c^2 x^2}} x^2}{40 c^3}-\frac {b \sqrt {1-\frac {1}{c^2 x^2}} x^4}{20 c}+\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {(3 b) \text {Subst}\left (\int \frac {1}{c^2-c^2 x^2} \, dx,x,\sqrt {1-\frac {1}{c^2 x^2}}\right )}{40 c^3} \\ & = -\frac {3 b \sqrt {1-\frac {1}{c^2 x^2}} x^2}{40 c^3}-\frac {b \sqrt {1-\frac {1}{c^2 x^2}} x^4}{20 c}+\frac {1}{5} x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {3 b \text {arctanh}\left (\sqrt {1-\frac {1}{c^2 x^2}}\right )}{40 c^5} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.06 (sec) , antiderivative size = 97, normalized size of antiderivative = 1.09 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {a x^5}{5}+b \sqrt {\frac {-1+c^2 x^2}{c^2 x^2}} \left (-\frac {3 x^2}{40 c^3}-\frac {x^4}{20 c}\right )+\frac {1}{5} b x^5 \sec ^{-1}(c x)-\frac {3 b \log \left (x \left (1+\sqrt {\frac {-1+c^2 x^2}{c^2 x^2}}\right )\right )}{40 c^5} \]

[In]

Integrate[x^4*(a + b*ArcSec[c*x]),x]

[Out]

(a*x^5)/5 + b*Sqrt[(-1 + c^2*x^2)/(c^2*x^2)]*((-3*x^2)/(40*c^3) - x^4/(20*c)) + (b*x^5*ArcSec[c*x])/5 - (3*b*L
og[x*(1 + Sqrt[(-1 + c^2*x^2)/(c^2*x^2)])])/(40*c^5)

Maple [A] (verified)

Time = 0.26 (sec) , antiderivative size = 141, normalized size of antiderivative = 1.58

method	result	size
parts	\(\frac {a \,x^{5}}{5}+\frac {x^{5} b \,\operatorname {arcsec}\left (c x \right )}{5}-\frac {b \left (c^{2} x^{2}-1\right ) x^{2}}{20 c^{3} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \left (c^{2} x^{2}-1\right )}{40 c^{5} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 c^{6} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}\)	\(141\)
derivativedivides	\(\frac {\frac {a \,c^{5} x^{5}}{5}+\frac {b \,c^{5} x^{5} \operatorname {arcsec}\left (c x \right )}{5}-\frac {b \left (c^{2} x^{2}-1\right ) c^{2} x^{2}}{20 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \left (c^{2} x^{2}-1\right )}{40 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, c x}}{c^{5}}\)	\(148\)
default	\(\frac {\frac {a \,c^{5} x^{5}}{5}+\frac {b \,c^{5} x^{5} \operatorname {arcsec}\left (c x \right )}{5}-\frac {b \left (c^{2} x^{2}-1\right ) c^{2} x^{2}}{20 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \left (c^{2} x^{2}-1\right )}{40 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, c x}}{c^{5}}\)	\(148\)

[In]

int(x^4*(a+b*arcsec(c*x)),x,method=_RETURNVERBOSE)

[Out]

1/5*a*x^5+1/5*x^5*b*arcsec(c*x)-1/20*b/c^3*(c^2*x^2-1)/((c^2*x^2-1)/c^2/x^2)^(1/2)*x^2-3/40*b/c^5*(c^2*x^2-1)/
((c^2*x^2-1)/c^2/x^2)^(1/2)-3/40*b/c^6*(c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)/x*ln(c*x+(c^2*x^2-1)^(1/2
))

Fricas [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 107, normalized size of antiderivative = 1.20 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {8 \, a c^{5} x^{5} + 16 \, b c^{5} \arctan \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) + 8 \, {\left (b c^{5} x^{5} - b c^{5}\right )} \operatorname {arcsec}\left (c x\right ) + 3 \, b \log \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) - {\left (2 \, b c^{3} x^{3} + 3 \, b c x\right )} \sqrt {c^{2} x^{2} - 1}}{40 \, c^{5}} \]

[In]

integrate(x^4*(a+b*arcsec(c*x)),x, algorithm="fricas")

[Out]

1/40*(8*a*c^5*x^5 + 16*b*c^5*arctan(-c*x + sqrt(c^2*x^2 - 1)) + 8*(b*c^5*x^5 - b*c^5)*arcsec(c*x) + 3*b*log(-c
*x + sqrt(c^2*x^2 - 1)) - (2*b*c^3*x^3 + 3*b*c*x)*sqrt(c^2*x^2 - 1))/c^5

Sympy [A] (verification not implemented)

Time = 3.42 (sec) , antiderivative size = 175, normalized size of antiderivative = 1.97 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {a x^{5}}{5} + \frac {b x^{5} \operatorname {asec}{\left (c x \right )}}{5} - \frac {b \left (\begin {cases} \frac {c x^{5}}{4 \sqrt {c^{2} x^{2} - 1}} + \frac {x^{3}}{8 c \sqrt {c^{2} x^{2} - 1}} - \frac {3 x}{8 c^{3} \sqrt {c^{2} x^{2} - 1}} + \frac {3 \operatorname {acosh}{\left (c x \right )}}{8 c^{4}} & \text {for}\: \left |{c^{2} x^{2}}\right | > 1 \\- \frac {i c x^{5}}{4 \sqrt {- c^{2} x^{2} + 1}} - \frac {i x^{3}}{8 c \sqrt {- c^{2} x^{2} + 1}} + \frac {3 i x}{8 c^{3} \sqrt {- c^{2} x^{2} + 1}} - \frac {3 i \operatorname {asin}{\left (c x \right )}}{8 c^{4}} & \text {otherwise} \end {cases}\right )}{5 c} \]

[In]

integrate(x**4*(a+b*asec(c*x)),x)

[Out]

a*x**5/5 + b*x**5*asec(c*x)/5 - b*Piecewise((c*x**5/(4*sqrt(c**2*x**2 - 1)) + x**3/(8*c*sqrt(c**2*x**2 - 1)) -
3*x/(8*c**3*sqrt(c**2*x**2 - 1)) + 3*acosh(c*x)/(8*c**4), Abs(c**2*x**2) > 1), (-I*c*x**5/(4*sqrt(-c**2*x**2
+ 1)) - I*x**3/(8*c*sqrt(-c**2*x**2 + 1)) + 3*I*x/(8*c**3*sqrt(-c**2*x**2 + 1)) - 3*I*asin(c*x)/(8*c**4), True
))/(5*c)

Maxima [A] (verification not implemented)

none

Time = 0.20 (sec) , antiderivative size = 131, normalized size of antiderivative = 1.47 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {1}{5} \, a x^{5} + \frac {1}{80} \, {\left (16 \, x^{5} \operatorname {arcsec}\left (c x\right ) + \frac {\frac {2 \, {\left (3 \, {\left (-\frac {1}{c^{2} x^{2}} + 1\right )}^{\frac {3}{2}} - 5 \, \sqrt {-\frac {1}{c^{2} x^{2}} + 1}\right )}}{c^{4} {\left (\frac {1}{c^{2} x^{2}} - 1\right )}^{2} + 2 \, c^{4} {\left (\frac {1}{c^{2} x^{2}} - 1\right )} + c^{4}} - \frac {3 \, \log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} + 1\right )}{c^{4}} + \frac {3 \, \log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} - 1\right )}{c^{4}}}{c}\right )} b \]

[In]

integrate(x^4*(a+b*arcsec(c*x)),x, algorithm="maxima")

[Out]

1/5*a*x^5 + 1/80*(16*x^5*arcsec(c*x) + (2*(3*(-1/(c^2*x^2) + 1)^(3/2) - 5*sqrt(-1/(c^2*x^2) + 1))/(c^4*(1/(c^2
*x^2) - 1)^2 + 2*c^4*(1/(c^2*x^2) - 1) + c^4) - 3*log(sqrt(-1/(c^2*x^2) + 1) + 1)/c^4 + 3*log(sqrt(-1/(c^2*x^2
) + 1) - 1)/c^4)/c)*b

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4828 vs. \(2 (75) = 150\).

Time = 0.79 (sec) , antiderivative size = 4828, normalized size of antiderivative = 54.25 \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\text {Too large to display} \]

[In]

integrate(x^4*(a+b*arcsec(c*x)),x, algorithm="giac")

[Out]

1/40*c*(8*b*arccos(1/(c*x))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*
x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c
^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 3*b*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/(c^6 + 5*c^6*(1/(c^2*x^2
) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)
^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) + 3*b*log(abs(sqrt(
-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/
(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6
*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) + 8*a/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x
^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) +
1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 40*b*(1/(c^2*x^2) - 1)*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c
^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x
) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)
^2) - 15*b*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1
/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^
6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 15*b*(1/(
c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 +
10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2)
- 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) - 10*b*sqrt(-1/(c^2*x^2) +
1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(
c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x)
+ 1)^10)*(1/(c*x) + 1)) - 40*a*(1/(c^2*x^2) - 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/
(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(
c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 80*b*(1/(c^2*x^2) - 1)^2*arccos(1/(
c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1
/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*
x) + 1)^10)*(1/(c*x) + 1)^4) - 30*b*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 +
5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) -
1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(
1/(c*x) + 1)^4) + 30*b*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^
2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x)
+ 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^
4) + 4*b*(-1/(c^2*x^2) + 1)^(3/2)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2
/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^
6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^3) + 80*a*(1/(c^2*x^2) - 1)^2/((c^6 + 5*c^6*(1/(c^2*x^2)
- 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^
6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 8
0*b*(1/(c^2*x^2) - 1)^3*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2)
- 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^
8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 30*b*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2
*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*
x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c
^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) + 30*b*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) -
1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 1
0*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^
5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 80*a*(1/(c^2*x^2) - 1)^3/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^
2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^
2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) + 40*b*(1/(c^2*x^2) - 1
)^4*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)
^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2)
- 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) - 15*b*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x)
+ 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(
1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c
*x) + 1)^10)*(1/(c*x) + 1)^8) + 15*b*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6
+ 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) -
1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*
(1/(c*x) + 1)^8) + 4*b*(1/(c^2*x^2) - 1)^3*sqrt(-1/(c^2*x^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1
)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*
x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^7) + 40*a*(1/(c^2*x^2) -
1)^4/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1
/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*
x) + 1)^10)*(1/(c*x) + 1)^8) - 8*b*(1/(c^2*x^2) - 1)^5*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x
) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/
(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) - 3*b*(1/(c^2*x
^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10
*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1
)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 3*b*(1/(c^2*x^2) - 1)^5*lo
g(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*
x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x)
+ 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 10*b*(1/(c^2*x^2) - 1)^4*sqrt(-1/(c^2*x
^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6
*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/
(c*x) + 1)^10)*(1/(c*x) + 1)^9) - 8*a*(1/(c^2*x^2) - 1)^5/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10
*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1
)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10))

Mupad [F(-1)]

Timed out. \[ \int x^4 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\int x^4\,\left (a+b\,\mathrm {acos}\left (\frac {1}{c\,x}\right )\right ) \,d x \]

[In]

int(x^4*(a + b*acos(1/(c*x))),x)

[Out]

int(x^4*(a + b*acos(1/(c*x))), x)